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We revealed a right-left asymmetry of the inter-well mean first passage times for the Brownian 
particles in a ratchet potential under internal white noise. We showed analytically and numerically 
that this asymmetry gives rise to the following phenomena: (i) nonstationary transport; (ii) local- 
ization or ultraslow (Sinai) diffusion for ratchets with a disorder in the asymmetry parameter of the 
potential; (iii) equilibrium "drift of labels" in circular ratchets. 
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It was demonstrated recently, that noncquilibrium fluc- 
tuations may support a stationary flux of the Brownian 
particles in a periodic potential with broken symmetry 
(ratchet). This fluctuation-induced transport has been 
of growing interest after Magnasco suggestion to apply 
a toy "ratchet and pawl" Feynman's model to riddle a 
" mystery" of the directed transport in cells Q . In sub- 
sequent studies , some related models were pro- 
posed to explain the nature of molecular motors , to 
develop a novel mass separation technique flffl , to design 
the Brownian engine with a negative resistance S, or to 
speculate how the ion channel junctions operate IjJ. Ac- 
cording to the second law, no stationary transport occurs 
under equilibrium fluctuations. To obtain a stationary 
motion, one has to break somehow the detailed balance. 
A possible way is to impose an additional colored noise 
[Q- 3^p|-|il| . One can also couple the system with nonequi- 
librium bath Jq], with external (e.g. chemical) degrees of 
freedom [^],|]||], or with a deterministic chaotic process 

In the present study we consider the motion of the 
Brownian particles in a ratchet potential (Fig. la) under 
an equilibrium white noise. We found a right-left asym- 
metry for the mean first passage times (MFPT) of the 
interwell transitions. This property of the MFPT entails 
the following phenomena: (i) nonstationary transport of 
the Brownian particles; (ii) ultraslow (Sinai) diffusion 
or "localization" in ratchets with a disorder (Fig. lb.); 
(iii) equilibrium "drift of labels" in circular ratchets. The 
directed nonstationary transport of the ensemble of the 
Brownian particles occurs as its peculiar route towards 
equilibrium Jl3| . If the process is interrupted at some 
stage (e.g. removing a particle) a net transport is ob- 
tained. This scenario of the nonstationary transport 
seems to be realistic for biological systems. 

Nothing is perfectly periodic in nature, and any peri- 
odic structure has imperfections and distortions. Some- 
times even weak disorder drastically changes a kinetic 
behavior of the system Jl4| ). If the ratchet potential 
is symmetric in average, but the asymmetry parameter, 
a = L-/Lo is stochastically distributed around a = 1/2 
(Fig.l), the mean displacement of the particles is zero. 



Surprisingly, the mean square displacement behaves as 
~ [log(i)] at large times t, contrary to ~ t dependence 
for the usual diffusion. At t 3> 1 [log(i)] jt — > 0, so the 
particle is, in fact, localized on the length scale of the 
ordinary diffusion. The "localization" may be effectively 
used in biological systems, or in technology, instead of 
trapping molecules in a deep potential well by forming a 
chemical bond. It costs much less energy to release the 
particle localized by the disordered ratchets, rather then 
waste energy, breaking the chemical bond. 

The directed nonstationary transport and ultraslow 
diffusion concern the Brownian motion on the infinite 
line. The other interesting phenomenon occurs in cir- 
cular ratchets. From one hand, the equilibrium flux in 
a ring is zero, from the other, the interwell MFPT for 
transitions to the right is not equal to the MFPT of the 
left transitions. As a result a "drift of labels" may be 
detected, provided that the particles may be labeled and 
their paths may be traced. The "drift of labels" is not 
a thermodynamical flux and and the detailed balance at 
equilibrium is preserved. 

To address these problems we consider the Brownian 
particle in a periodic ratchet potential (Fig. la), subjected 
to the white noise, £_{t), < >= 2DS(t-t'). Here 

D = kT/j is the diffusion, and 7 is the friction coeffi- 
cient. The Langevin equation for the overdamped case 
reads: 



(i) 



with f(x) = —7 1 dU/dx, and U{x) being the ratchet 
potential. The corresponding Fokker-Plank equation is 



d t p = -d x (fP) + do xx p 



(2) 



The Brownian motion in the piecewise linear potential 
may be effectively studied in terms of the MFPT [ p3| |. 
The escape time from a starting point, xq, in the well is 
given as a solution of the equation p6| : 
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The final point is either L + , if the particle escapes to the 
right, or — if it overcomes the barrier in the oppo- 
site direction. The averaging of the escape time over Xq 
gives T+ and T_, the MFPT for the right and for the 
left transitions respectively. To find T+, we use the con- 
tinuity of T'(xo), and boundary conditions T(L + ) = 0, 
T'(-L-) = 0, which are self-consistent at t±/T± <C 1. 
Here t± = r yL±/kTa± are the intrawell relaxation times, 
and a± = Uo/kTL± characterize the force on the slopes. 
The boundary conditions are reversed to obtain T_ . We 
consider the case when exp(Uo/kT) 3> 1 and t±/T± -C 1, 
i.e., the particle has enough time to thermalize inside the 
well before it jumps over the barrier. Straightforward 
analysis shows that the solution of Eq.(||) is a sum of 
few terms fll^] . The leading term is independent on xo, 
and is of the order of exp(Uo/kT) 3> 1. The other terms 
|jl7f may be neglected everywhere exept the regions in 
the very vicinity of the tops. If the distribution inside 
the well is sharply peaked around the bottom, these re- 
gions do not noticeably contribute to the MFPT. Par- 
ticularly, this is true for the equilibrium distribution if 
U /kT > 0.1. Finally, we obtain 



To describe the nonstationary transport we introduce the 



T± = a ± 1 e u "/ kT /D 



«o 



(4) 



where aZ = a+ + aZ ■ This equation shows that 
T + T_ for the potential with asymmetry. It is worth 
noting that the left-right asymmetry of the MFPT for the 
ratchet under the white noise occurs as the preexponen- 
tial factor. For the colored noise this factor is overshad- 
owed by the asymmetry that appears in the exponent 
. The left-right asymmetry of the MFPT described 
by Eq.(Q) was confirmed in our numerical studies per- 
formed for the equilibrium intrawell distribution. To get 
the required white noise process (which is effectively col- 
ored due to the finite time step) we decreased the time 
step until the results became insensitive to it further de- 
crease in a range of two decades (the time step varied 
from 2 10 -4 to 10 -6 , and was roughly two orders of mag- 
nitude smaller, then that, used in p8|). We found that 
the simplest Euler method (see e.g. [18|) is quite effec- 
tive for the piecewise linear potential. 

As a consequence of the left-right asymmetry in the 
MFPT, the directed motion appears before the system 
reaches the stationary probability distribution, P eq (x) oc 
exp(— U(x)/kT), which is a current free, if no density 
gradients present in the system. Thus, any initially lo- 
calized distribution, say P(x, t\ xq, to) = S(x — Xo)(t — to), 
evolves to the P eq (x), defined on the ratchet, with left 
and right boundary — L^, L m L$. In the sequence 
of the time scales, that characterize the Brownian mo- 
tion in the overdumped case, the shortest time scale cor- 
responds to the intrawell relaxation time t±. The next 
time scale is related to the interwells transition times, T±. 
Finally, ~ (Loo/ Lq)T± defines the global relaxation 
time. Under condition exp(Uo/kT) ^> 1, the following 
time-scales hierarchy holds, t± <C T± <C T^, and the 
nonstationary transport takes place for T± <C t <C Too- 



population of the k-th well as Ufc(t) = f Lk 
and write the kinetic equation: 



LoO+l) 



P{x, t)dx, 



W^n k+1 + W+n k -i - (W- + W+) n k 



(5) 



where W± give transition rates, which may be estimated 
as W± ~ TZ 1 fl9|] . Performing transformation from n k 
in the discrete space to n(x) in continuous space coordi- 
nate, we obtain the coarse-grained equation: 

d t n(x, t) ~ -Vd x n(x, t) + Dd xx n(x, t) (6) 

with V = LoiT- 1 - TZ 1 ) and D = \LKTl 1 + TZ 1 )- 
More rigorous derivation of the coarse-grained equation, 
based on the separation of different time-scales |20| yields 
the same result for t ^> T±. 

Eq.(||) shows that the ensemble of Brownian parti- 
cles moves with a constant velocity V, which increases 
with enhancing asymmetry of the potential as V = 
(U Q /kT) 2 e-^ kT HD/L )(j^ - i) (Fig.2). The in- 
crease of the distribution width with time is described 
by the renormalized diffusion coefficient D. 

We apply the coarse-grained description for the case 
of the external field U ex (x) = —Ex. Calculations of the 
MFPT for the small fields, EL /U Q < 1, yield Eq.(|), 
with a± substituted by a± ± E/kT, and Uo/kT substi- 
tuted by (Uq ± EL±)/kT. In the linear approximation 
for E, the average velocity reads: 



V E = V + (D/kT)E 



(7) 



where V and D are the field-free values. The top insert in 
Fig.2 depicts dependence of the average velocity upon the 
external force. The particles move uphill, if the field does 
not exceed the critical value, E c — (2kT/L )(l — 2a). 

The kinetic behavior of the system in the vicinity of 
the critical field, i.e. near Ve = drastically changes 
under (even weak) disorder. The disorder may be ei- 
ther in the values of E, or in the values of T±, i.e. 
in the asymmetry parameter a. For such a "quenched 
disorder" Q, one has stochastically distributed tran- 
sition probability rates, W + ~ TZ 1 and VF_ ~ TZ 1 
(probabilities for the right and for the left "jumps"). If 
(log(W + )) — (log(VF_)), the average displacement of the 
particle is zero, while the average square displacement 
behaves asymptotically as ~ (logi) 4 [ pT|j2l| . This be- 
havior is called as ultraslow (Sinai) diffusion, and may be 
regarded as a localization in the length scale of the ordi- 
nary diffusion Jl4|]. If the external field is zero, E = 0, 
the ultraslow diffusion occurs in ratchets, which are sym- 
metric in average. 

To prove the existence of the ultraslow (Sinai) diffu- 
sion, we applied an approach commonly used to analyze 
this phenomenon pl[ . Namely, we consider a motion 
of a particle in a segment of N = 20 ratchets with a 
distributed uniformly between 0.2 and 0.8, so that the 
average length of the left slope was equal to the average 
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length of the right slope. The reflecting left boundary and 
adsorbing right boundary were used. We calculated nu- 
merically the distribution P(Tjv) of times TV, when the 
particle, initially located at the left boundary, reached 
the adsorbing boundary. For the ordinary Id diffusion, 
the Tv S> 1 asymptotics reads P(Tjv) ~ T N 3 ^ 2 , while 
for the Sinai diffusion P(T N ) ~ T^ 1 The numer- 

ical results shown in Fig. 3 illustrate that the ultraslow 
diffusion does take place in ratchets with the disorder. 

The transition probability rates W±, as well as the 
kinetic equation ( ^) refer to the time coarse-grained de- 
scription and may be unambiguously applied only for 
<>T±. To analyze the detailed balance or the thermo- 
dynamic flux J(x) = limAt^o(^+ ~ N-)/At, one needs 
more refined approach, since At <C T±. Here N + and 
denotes correspondingly the number of particles that 
have crossed the plane located at x in positive and nega- 
tive directions during the time interval At. To show that 
the asymmetry of the T± (x) does not imply the non-zero 
equilibrium flux, we write for the flux, e.g. at x = L + (see 

Fig.l): J(x,t) = (l/At)[/f+_ P(x,t) 6{At-T+{x))dx- 

f£ ++L ° P(x,t) 6{At - T..(x))dx\, where the unit step 

function 9(t) guarantees that particles have crossed the 
pla ne within the At. Since At —> 0, we expand T±(x) 
|7| near x = L+: T ± {x) =| L+ - x \ e u °/ kT /a D, As a 
result, a zero flux is obtained for the equilibrium proba- 
bility density P eq {x). 

In spite of a lack of the thermodynamic current at 
equilibrium, very unusual equilibrium fluctuations may 
be observed in circular ratchets. If particles may be la- 
beled and if their paths may be traced, one observes that 
the average number of times that particle crosses some 
cross-section clockwise, (5+), is not equal at equilibrium 
to the average number of the counter-clockwise cross- 
ings, (S-). Therefore, clockwise (or counter-clockwise) 
"drift of labels" takes place in equilibrium. Obviously, 
this is not a thermodynamical current, and no useful 
work may be extracted. We quantify the "drift of labels" 
by S — ({S+)/{S-) — 1). We also introduce the mean 
times of the clockwise turn over the ring, Tc+, and of the 
counter-clockwise turn, Tc—. If, again, exp(Uo/kT) 3> 1 
and t±/T± <C 1, a straightforward solution of Eq.(||) 
gives for the ring of M ratchets: 

Tc± = ( M(M-l) + M\ 

Dciq \ 2ao a± J 

and we estimate S as (T^j/T^ - 1) = (M 4 + "^ 2et ■ The 
"flux of labels" disappears at large M, and thus it is 
an intrinsically microscopic phenomenon. Our numerical 
studies of the Brownian motion on the ring, of M = 4 
ratchets with a = 0.8 confirmed that the average current 
ceased to zero, while S relaxed to 0.29, which is in a good 
agreement with the theoretical estimate S = 0.35. 

In conclusion, we studied the Brownian motion of par- 
ticles in a ratchet potential under the internal white 
noise. We found that the inter-well mean-first passage 



times have a right-left asymmetry for the equilibrium 
and near-equilibrium intrawell distributions of particles. 
We have shown analytically and numerically that this 
asymmetry gives rise to the following effects: (i) nonsta- 
tionary transport; (ii) localization or ultraslow (Sinai) 
diffusion under disorder in the parameter of asymmetry 
of the ratchet potential; (iii) unusual equilibrium fluctua- 
tions in circular ratchets, which may be called as a "drift 
of labels" . 

We are indebted to M. Bier, R. Kapral, and A. Male- 
vanets for valuable comments and discussions. 
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Fig. 1 Ratchet potential, L = L + + L- (a); ratchet 
with a disorder of the asymmetry parameter (b). 

Fig. 2 The average displacement versus time for the 
ensemble, initially localized at x = 0. Bottom insert: the 



average velocity, V, versus the asymmetry parameter, a. 
Top insert: V versus the external force, E. The dots 
- simulations, the lines - theory; kT — 0.085, a = 0.8, 
7=1, I,o = l, U = l. 

Fig. 3 The distribution of escape times for the seg- 
ment of 20 ratchets with reflecting left boundary. Bottom 
curve - the symmetric ratchet, a = 1/2. Upper curve - 
the ratchet with a disorder; a flat distribution of a was 
used. 
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